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DAMPING— IN-ROLL  CALCULATIONS  FOR  SLENDER  SWEFT-^BACK 
WINGS  AND  SLENDER  WING-^ODY  COMBINATIONS 
By  Harvard  Lomax  and  Max.  A.  Heaslet 


SUMMARY 

The  damping— in-roll  parameter  C^p  is  calculated  theoretically 
for  triangular  wings  on  cylindrical,  bodies  and  for  a  class  of  wings 
with  swept-back  plan  forms.  The  analysis  is  based  on  the  usual  assump¬ 
tions  of  linearized  compressible-flow  theory  together  with  the  added 
restrictions  that  at  the  free— stream  Mach  number  Mq  the  product  of 
1_Mq2  and  the  streamwise  velocity  gradient  is  small.  The  accuracy  of 
the  results  is  tested  by  a  comparison  with  the  exact  solution  of  the 
linearized  equation  for  a  triangular  wing  and  for  such  plan  forms  is 

shown  to  increase  as  reduced  aspect  ratio  (v|l-M02|  times  aspect 
ratio)  decreases. 


*  INTRODUCTION 

Theoretical  advancements  in  the  study  of  load  distributions  over 
three-dimensional  wings  in  compressible— flow  fields  have  been  achieved 
almost  entirely  under  assumptions  leading  to  the  linearization  of  the 
basic  differential  equations  of  flow.  In  this  manner  the  analysis  is 
resolved  in  general  into  the  problem  of  determining  solutions  of  elliptic 
and  hyperbolic  partial  differential  equations  in  three  dimensions,  the 
character  of  the  equations  being  fixed  by  the  magnitude  of  the  free-stream 
Mach  number  Mq*  It  is  possible,  however,  to  study  also  a  class  of  prob¬ 
lems  associated  with  the  parabolic  form  of  the  potential  equation;  namely, 
cases  for  which  the  product  of  1-M0S  and  the  chordwise  perturbation- 
velocity  gradient  is  small  in  comparison  with  the  gradients  in  their 
respective  directions  of  the  perturbation  velocities  perpendicular  to  the 
chord.  Such  a  method  was  developed  in  reference  1  for  small  values  of 
aspect  ratio  and  was  used  in  reference  2  for  swept-back  lifting  surfaces 
at  sonic  speeds.  In  the  present  report  the  method  will  be  referred  to 
as  a  slender-wing  theory  wherein  the  term  slender  infers  that  the  ratio 
of  the  reduced1  span  to  the  over— all  length  is  small.  Hence  the  results 

deduced  span  is  defined  as  ft  times  the  span,  where  ft  is  equal  to 
fj  |  1-Mq2  |  .  In  general,  the  term  reduced  implies  that  any  parameter 
which  it  modifies  is  multiplied  by  the  factor  ft.  • 
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presented  will  apply  most  closely  when  the  wing  is  long  in  terms  of  its 
span,  or  when  the  free— stream  Mach  number  is  close  to  unity  and  the  wing 
and  hody  do  not  violate  the  assumptions  made • 

Previous  investigations  (references  3,  5,  and  6)  corresponding 

to  this  method  of  approach  include,  respectively,  a  complete  analysis  of 
all  the  stability  derivatives  of  a  low-aspect-ratio  triangular  wing,  a 
lifting  triangular  wing  with  an  arbitrary  body  of  revolution,  a  lifting 
triangular  cruciform  combination  on  an  arbitrary  body  of  revolution,  and 
a  lifting  swept-back  constant-chord  wing.  In  order  to  assess  the  accu-  • 
racy  of  the  results  presented  in  these  reports,  it  is  necessary  to  com¬ 
pare  their  approximate  solutions  with  the  exact  solutions  of  the  linear¬ 
ized  equation.  For  the  case  of  the  rolling  triangular  wing,  this  can 
easily  be  done  since  the  exact  solution  has  been  derived  in  reference  7* 
The  results  of  such  a  comparison  are  presented  in  figure  1  which  indi¬ 
cates  that  in  this  case  the  value  of  (BCjp  for  the  approximate  theory 
can  be  useful  up  to  reduced  aspect  ratios  as  high  as  three. 

The  object  of  the  present  report  is  three— fold:  First,  to  find  the 
effect  on  the  damping— in— roll  parameter  C^p  of  adding  a  body  of  revolu¬ 
tion  to  a  triangular  wing  in  order  to  extend  the  knowledge  of  wing-body 
interference  into  the  field  of  lateral  stability;  second,  to  find  CZp 
for  a  particular  swept-back  wing  plan  form;  and  third,  to  show  how  the 
damping  due  to  roll  can  be  found  for  swept  wings  with  arbitrary  trailing 
edges  or  how  trailing  edges  can  be  calculated  from  prescribed  span  load 
distributions.  The  examples  to  be  given  provide  sufficient  details  to 
indicate  how  other  cases  can  be  calculated  numerically. 

A  list  of  important  symbols  is  given  in  Appendix  A. 

RESUME  OF  THE  METHOD 

Under  the  assumption  that  (l-M02)$^  is  small  as  compared  to  Oyy 
and  the  equation  for  compressible  fluid  motion,  either  subsonic 

or  supersonic,  becomes 


® IY  +  ®ZZ  “  0  (la) 

where  $  is  the  perturbation-velocity  potential  and  X  is  the  free— 
stream  flow  direction.  The  resulting  differential  equation  (la)  shows 
that  although  $  can  be  a  function  of  X  its  value  is  a  consequence  of 
boundary  conditions  given  along  lateral  strips.  We  will  seek  solutions 
to  equation  (la)  for  which  first,  the  perturbation  velocities  vanish  at 
infinity;  second,  there  is  no  discontinuity  in  potential  except  across 
a  lifting  surface  or  its  trailing  vortex  sheet;  and  third,  there  is  no 
discontinuity  in  the  vertical  velocity  anywhere  in  the  field  (i.e.,  no 
airfoil  with  thickness).  It  is  possible,  by  applying  the  two-dimensional 
form  of  Green fs  theorem,  to  write  the  solution  to  equation  (la)  which 
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satisfies  all  of  the  aforementioned  conditions  and  relates 

(Oz  in  the  Z  =  0  plane)  to  Avt0  (the  jump  in  sidewash,  G>y  in  the  Z  =  0 

plane).  Under  the  added  restriction  that  no  point  exists  in  or  on  the 

flow  field  such  that  ^  0,  where  €  is  the  radial  distance  from 

— Be 

the  point,  this  solution  can  he  written2 


Wto  2n 


,B(X)  Avt  dYi 


-B(X) 


The  integral  equation  (2a)  will  be  solved  for  the  case  of  rolling 
wings.  First,  however,  it  is  convenient  to  introduce  nondimens ional 
lengths  so  that  the  results  can  be  generalized  as  far  as  possible.  A 
satisfactory  nondimens ional  form  is  obtained  by  using  the  equations 


C  tan  0 


=  C  v-fc  tan  6 
=  C  w-t-  tan  9 


=  C  v^  tan  9 


C  tan  0 


=  p  C2  tan2  9 


where  9  is  the  semivertex  angle  of  the  wing,  p  is  the  rate  of  roll  in 
radians  per  second,  C  is  the  root  chord,  u^,  v^.,  w^, and  vr  are  true  of 

velocity  components,  and  u,  v,  w,  and  vr  are  the  transformed  components. 
In  addition  cp(x,y,z)  replaces  $(X,Y,Z). 


TRIANGULAR  WING  WITH  BODY 


Consider  a  wing  with  an  unswept  trailing  edge  mounted  on  a  body  of 
revolution  as  shown  in  figure  2.  Applying  the  transformations  given  by 
equations  (3)  to  equation  (la),  we  find  that  in  every  plane  perpendicular 
to  the  x  axis  it  is  necessary  to  satisfy  the  equation 


2This  equation  is  developed  in  more  detail  in  reference  2. 
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subject  to  the  boundary  conditions 


it 

0 

for 

r  =  a,  0  < 

t  <  2* 

w0  =  ny 

for 

z  =  0,  a  < 

|y|  <b 

y  =  y  =  0  for  r=a>,  0<ty<2lt  J 

where  a(x)  and  b(x)  are  replaced  by  a  and  t  for  convenience  in 
notation  and  where  r,t  are  polar  coordinates  in  the  yz  plane  and  p 
is  related  to  the  rate  of  roll  by  equation  (3).  These  boundary  condi¬ 
tions  are  not  of  a  type  such  that  they  can  be  substituted  directly  into 
equation  (2a),  since  one  condition  specifies  a  radial  component  which 
is,  moreover,  out  of  the  z  =  0  plane.  However,  equation  (1)  is  simply 
Laplace f s  equation  so  that  by  introducing  complex  variables  it  is  not 
difficult  to  transform  the  wing  and  body  section  (fig.  3(a))  to  a  strip 
along  the  y  axis  (fig.  3(b)).  Such  a  mapping  function  is  provided  by 
the  Joukowski  transformation,3  which,  if  the  subscript  1  denotes 
conditions  in  the  transformed  plane,  can  be  written: 


(5) 


where  |  and  are  given  by  the  expressions 

i  =  y  +  iz 

.  ► 

§1  =  yi  +  izi 


(6) 


By  means  of  equation  (5)  the  wing  in  the  |  plane  (the  portion  of 
the  axis  for  which  a<|y|  <b,  fig.  3)  transforms  into  the  section  o| 
the  yi  axis  in  the  |i  plane  for  which  2a  =  ax  <|  yL  |  <bi  -  3  +  b  * 
Further,  the  body  in  the  |  plane  (the  curve  satisfying  the  equation  r=a) 
transforms  into  the  part  of  the  71  axis  for  which  0  ^|yi|  <ax  -  <>a. 

In  this  manner  the  boundary  on  which  the  data  are  specified  has  been 
transformed  so  that  it  lies  entirely  along  the  yx  axis.  It  is  next 
necessary  to  inspect  what  the  boundary  conditions  are  m  the  61  plane. 

From  the  basic  theory  underlying  the  use  of  complex  variables,  induced 
velocities  in  the  physical  and  transformed  planes  are  related  by  the 
expression 

v  —  iw  =  (vx  —  iWi)-i-' 

d| 

3See,for  example,  "The  Elements  of  Aerofoil  and  Airscrew  Theory," 

Cambridge  University  Pres6,  19^3,  "by  H.  Glauert. 


(7) 


t* 


I 


9 


» 
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from  which.,  since  in  polar  coordinates 


cL|i  = 
dl 


(0 


cos  2  i|r 


+  i 


(0 


sin  2 


(8) 


it  follows  Immediately  that 

2 


V  =  v± 


vr  = 


2  i i 

+  “i  sin 

2  \|r 

-  (fj  sin 

sin 

1 - 1 

C\J^ - 

1 

1  H  1 

1 - 1 

-> 

>  (9) 


+  2  Wi 


0 


sin  \|r 


Through  the  use  of  equation  (9),  the  boundary  conditions  given  by 
equation  (4)  for  the  |  plane  can  be  transformed  to  the  plane  with 

the  result  (since  \|r  =  0  or  n) 


*1 

=  0 

for 

0  <  |yi[< 

Wt 

=  £  • 

yi 

2yi2-ai2 

4 

l7M 

V  yi2-aa2 

Vl 

= 

=  0 

for  tx  = 1 

The  problem  has 

now  became  one 

y  do) 


ary  conditions  given  by  equation  (10).  But  this  problem  can  be  immediately 
solved, provided  we  can  invert  the  integral  equation  (2a),  which  in  the  new 
notation,  becomes 


vi(yi) 


_  i.  f 'bl  Av1(yg)dyg 
2*  tLt)1  yi-ys 


(2b) 


Note  that  the  subscript  o  (indicating  conditions  in  the  Zx  =  0  plane) 
has  been  dropped  when  in  the  |  x  plane*  This  avoids  a  cumbersome 
notation  and  should  cause  no  confusion,  since  in  the  plane  only 

conditions  on  the  jx  axis  are  to  be  considered. 
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Under  the  condition 


'x  wx(yx)d 7i 


,  2  _  2 
'a  — yi 


reference  8  gives  an  inversion  to  equation  (2b)  which  can  be  written 


Wi(y2)dy2  cQ 

'  ■  '  '»'■«  " ,IW 

( yi— y2  )*J  "bi2— y22  */  bx2— yx‘ 


where  cQ  is  an  arbitrary  constant. 

An  inspection  of  vx  as  given  by  the  boundary  conditions  in  equation 
(10)  shows  that  it  is  an  odd  function  of  yx  and  hence  will  satisfy  equa¬ 
tion  (11).  Placing  wx  in  equation  (12)  and  integrating  gives 

Avi(yi)  =  _  J b^-yx2  (|  +  arc  cos  - 

2  v  ax2— yx2 


VYi  ln  yx»/bx2-ax2  +  ax»/bxg-yxg  + 
rt  yx  */b x 2— ax2  -  ax  */b12-y12  V/bI2^yj? 


2  „  2 


where  c  =  -1  for  0  <|yx|<ax,  and  cr  =  +1  for  ax  <|yx  |<  bx. 

Now  to  find  Acpx  on  the  zx  =  0  plane,  the  relation 

pY  1 

Acpi  =  /  Avx(y2)<ly2  (1^) 

-bx 

in  which  Avx(y2)  as  given  by  equation  (13)  is  used  with  the  result  that, 
for  ax<  yx<bx. 


|jl  /_  2  ax 

"5  f  1  +  5f  ar°  cos 


2  2 
-yi 


(yi2-ai2)  arc  cosh gk  /A^| 

hx  Vyx^-ai2 


¥ 
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where  the  constant  cQ  has  been  evaluated  by  the  conditions  that  Acpx 
must  vanish  when  yx  =  bx.  This  can  be  written  in  terms  of  the  original 
coordinates  (fig.  2(b))  and  boundary  conditions (equation  (4) ) ,  by 
retransforming  in  accordance  with  equation  (5).  Thus,  for  a<y<b. 


Acp0 


cos 


2ab 


b2 


+  a‘ 


)(-#)/ 


a 

b2 


-y2- 


a* 

y2 


arc 


cosh 


(y^-fa2)  (h2~a2) 
(y2 —a2)  (b2+a2) 


(16) 


V 


Equation  (16)  represents  the  solution,  in  terms  of  the  velocity  poten¬ 
tial,  for  a  rolling  wing  with  an  unswept  trailing  edge  on  a  body  of  revolu¬ 
tion]  both  a  and  b  are  functions  of  x.  The  loading  coefficient  ^ 
can  be  determined  by  taking  the  partial  derivative  of  equation  (16) 
with  respect  to  x  in  order  to  find  Au,  and  then  using  the  relation 


Ap  _  2Au£  _  2Au 


CV. 


(17) 


However,  it  is  not  necessary  to  find  the  loading  over  the  wing  if  interest 
is  limited  to  the  determination  of  the  total  rolling  moment  L.  If  a  is 
independent  of  x  (i.e.,  the  body  is  a  circular  cylinder),  then  r  =  (Acp)^ 
and  the  equation 


E. 


dZ  =  p  V0  T  dY 

where  dZ  is  the  increment  in  span  loading,  can  be  used  to  express 

dl  =  YdZ  =  p  V0  Y  T  dY 


Thus 


or 


L  =  2p  Vn  y^Y(A$0)m  -c  dY 


T.E. 


n® 

L  =  2p  V0  C2  tan2  6  /  y(Atp0)_  _  <4/ 

J  T.E. 


(18) 


where  _(A%)t>e  is  the  value  of  A%  at  the  trailing  edge  or  where 

can  deiiJlgUS  *  Su*?titutin6  equation  (16)  in  equation  (18),  one 
can  derive,  after  some  rather  involved  integration,  the  result 
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=  _  cf  tan4  6  E 
PV0 


9{h 


(a2+s2) 


- - axe  tan  — 

2  a 


JL  (s2-a2)  ^a^-6a2s2+s4  ^  arc  tan  ^  — 


&£  +  I  4  ( s2— a2 ) 


o,2  1 


2  s 


i 


(19) 


The  damp ing-in-r oil  coefficient  Cjp  based  on  the  wing  area  including 
the  part  of  the  wing  masked  by  the  body  can  now  be  written,  if  both 
sides  of  the  equation  are  multiplied  by  P  =  v  1 1-M02 !  >  in  "the  form 

_2 

+  ZR  ( 1-R2 )  (E  -6R2+l)  arc  tan  ^  - 
«2r4  +  E2(i^2)2  |  (20) 

where  E  =  -  =  -  (i.e.,  the  ratio  of  the  body  diameter  to  the  wing  span) 
s  3 

and  Ar  is  reduced  aspect  ratio  which  is  also  based  on  the  total  wing 
area  including  that  part  masked  by  the  body. 

When  E  equals  zero,  equation  (20)  represents  the  damping  in  roll 
for  a  wing  without  body.  This  value,  which  will  be  designated  (Czp)w, 
is  given  by 


•>%  ’  -  &  { 


(1+E  )  arc  tan  ± 


P 


(Cip)v  - 


jcAj. 

32" 


(21) 


and  agrees  with  the  value  given  in  reference  3-  Tbe  ratio  of  equation 
(20)  to  equation  (21)  represents,  finally,  the  effect  of  adding  a  circu- 
lar  cylinder  about  the  line  of  symmetry  of  any  pointed,  low-reduced-aspect- 
ratio  wing  with  an  unswept  trailing  edge.  This  effect  is  shown  in  figure 
The  figure  shows  that  the  damping  in  roll  is  increased  a  maximum  of  about 
4  percent  at  a  value  of  E  equal  to  0.28. 


SWEFT-BACK  WING 

The  study  of  the  swept-back  wing  will  be  divided  into  two  parts. 

The  first  of  these  sections  will  be  concerned  with  the  development  of  a 
special  plan  form  which  is  relatively  easy  to  treat  analytically,  while  the 
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second  section  will  be  devoted  to  the  examination  of  plan  forms  with 
arbitrary  trailing  edges  or  with  arbitrary  span  load  distributions. 


Special  Plan  Form 

Consider  a  wing  with  a  plan  form  similar  to  that  shown  in  figure  5. 
The  analysis  of  such  a  wing  is  best  presented  by  considering  separately 
the  three  regions  shown  in  figure  5(t>)  which  correspond, respectively, to 
values  of  x  for  which  x<l,  l<x<s,  and  s£x. 

Begion  I  (x< l).—  Since  the  loading  in  each  region  is  independent  of 
the  flow  at  all  points  downstream  of  it,  the  solution  for  the  pressures 
in  region  I  is  the  same  as  that  for  a  triangular  wing.  Hence  region  I 
is  simply  a  special  case  of  the  preceding  section  for  which  the  body 
diameter  B  equals  zero. 

Begion  II  (l£x£  s).—  Although  each  region  is  certainly  independent 
of  subsequent  ones,  the  vorticity  of  upstream  regions  must  be  transported 
downstream  through  all  subsequent  transverse  planes.  This  increases 
somewhat  the  complexity  of  the  problem. 


In  order  to  find  the  solution  for  the  loading  in  region  II,  it  is 
again  necessary  to  satisfy  equation  (l),  but  this  time  subject  to  the 
boundary  conditions: 


where5  a(x)  and 
and  b. 


w0  =  ny  for  a  <  |y  |<  b 
u  =  w  =  0  for  r=oo,0<i|r<2Jt  V 
Auq  =  0  for  0  <  |y  |<a 

✓ 

b(x)  are,  as  in  the  preceding  section,  replaced  by 


(22) 


a 


It  is  apparent  that  no  transformation  of  the  plane  is  necessary  for 
the  use  of  the  integral  equation  (2),  since  the  boundary  conditions  are 
already  specified  in  the  z  =  0  plane.  Making  use  of  the  fact  that  A<p 
must  be  an  odd  function  of  y,  and  hence  its  derivative  Av0  even, 
equation  (2b)  can  be  written 


b2 

Zo  =  u  =  _  r  Q(x^i)dni 

y  2k  J 


(23) 


5Whereas  in  the  wing-body  section  a(x)  represented  the  radius  of  the 
body,  here  a(x)  represents  the  distance  from  the  x  axis  to  the 
trailing  edge.  Since  the  two  sections  are  entirely  separate,  this 
should  cause  no  confusion. 
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where 


and 


*\ 


G(x,tj) 


AvQ(x,y) 

y 


> 


(24) 


Before  the  solution  to  equation  (23)  can  be  obtained,'  however,  the  behav 
ior  of  the  flow  in  the  vortex  wake,  the  region  for  which  -a<  y<  a,  must 
be  discussed. 


Since  an  increase  in  Acp0  corresponds  to  a  loading  A%  at  any 
point  behind  the  wing  must  be  the  same  as  Acp0  at  the  trailing  edge  for 
the  same  y.  If  the  equation  of  the  trailing  edge  is  written  in  either 
of  the  forms 


y  =  a(x) 
x  =  a*(y) 


(25) 


then  Acp0(x,y)  equals  Ap0(a*,y)  for  x>a*.  Taking  the  partial  derivative 
of  Acpo(a.*>y)  with  respect  to  y  gives 


dAcp0(a*,y)  =  &^Pp(a*,y)  da* 

c)y  c)a*  Sy 


dA%(x,y) 

dy 


x=a* 


which  can  also  be  written  in  the  form 


dAqb(a*,y)  =  V:  AP (»**?)  +  Avc(a*,y)  (2 6) 

dy  2  dy  q 


Applying  the  Kutta  condition,  the  value  of  the  load  coefficient  at  the 
trailing  edge  4E.(a*,y)  must  be  zero.  Equation  (23)  can  now  be  written 

1  pa2G*(ni)d-ti  _  2_  rb22iljM}^L  (27) 

2rt  J  T|— vji  2k  J  ti-tix 

o  a 

where  G*(t]x)  =  G(a*,tu). 


The  value  of  G*(t]i)  depends,  of  course,  on  the  equation  of  the 
trailing  edge.  It  is  convenient  at  this  stage  of  the  analysis  to  choose 
the  value  of  Acp0  in  the  wake  (since  applying  the  Kutta  condition  to 


equation  (26)  yields 


G*(m)  =  J 


^A<p0(x,y)  " 
-  dy 


x=a*  ’ 


giving  Acp0  is  equivalent 


* 


* 


MCA  TN  1950 


11 


to  specifying  G*(t|i) )  and  find  the  resulting  trailing  edge  rather  than 
to  proceed  in  the  more  straightforward  way  of  first  completely  specifying 
the  wing  plan  form.  However,  if  Acp^  is  expressed  in  the  form 

— ^  =  Bxy  +  Bay2  +  •  •  •  +  Bmym  (28a) 

Hi 

it  will  be  shown  that  both  the  trailing  edge  and  the  final  value  of 
can  be  given  in  terms  of  the  Bn's,  and,  therefore,  quite  general  types 
of  plan  forms  can  be  analyzed.  Our  purpose  is  to  present  in  some  detail 
the  simplest  example  in  this  section  and  to  develop  the  general  discussion 
later . 


Figure  6  shows,  the  variation  of  ATo  in  the  wake  behind  a  triangu¬ 
lar  wing  (equation  (15)  with  a  =  0,  b  =  l) .  A  reasonable  first  choice  of 
Acp0  in  the  wake  of  a  swept— back  wing  would  seem  to  be  the  straight-line 
continuation  of  the  slope  at  the  origin  of  figure  6.  This  leads  to  the 
equations 

Bj.  =  1,  B2  =  B3  =  •  .  .  =  0 

so  that  equation  (28a)  becomes 

G*(iu)  =  —  M-/  */tJT  '  (28b) 

6 

If  G(x,t| x )  is  defined  by  the  relation 


G(x,tii) 


+ 


(29) 


where  H(x)  is  a  function  to  be  determined  independent  of  t|i> 
equation  (27)  can  be  written 


M- 


2rt 


dq 


(n^liV  Tli 


H(x) 

2 


1_  r°  G1(x,T)i)dq1 

2rt  J& 2  n-Tii 


then 


(30) 


It  is  now  possible  to  invert  equation  (30)  means  of  equation  (12), 
wherein  a2  now  replaces  — b^  and  b2  replaces  b^.  First  we  notice 
that  the  condition  imposed  by  equation  (ll)  is  satisfied  if 


0 

The  choice  of  such  a  definition  is,  of  course,  not  obvious.  Its  advis¬ 
ability  depends  entirely  on  the  simplification  which  it  introduces 
into  the  subsequent  analysis. 


12 
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foil 

(T)-ni)-/ni 


H(x) 

2 


*) 

*/  (b2-T))(Tj-a2) 


=  0 


or,  upon  integrating,  that 

EW  =%.(»_£) 

where  K  is  the  complete  elliptic  integral  of  the  first  kind  of  modulus 
k  =  2:.  The  solution  for  Gx(x,t))  in  the  interval  a  <|y|<b  "becomes 


Gi(x,t|)  =  ^  «/(b2-r])(iy 


.b2 


1 


iL 

2rt 


dri2 


(T11-T12)  VruT 


*1 


( T]— 11 1 )  v  (t  -ill) (rji-a  ) 


The  single  integral  evaluates  to  give  zero,  and  the  order  of  the  double 
integral  can  he  reversed  and  integrated  with  respect  to  r\x  so  that 
finally  there  results 


Gi(x,n)  ("b^HTl-a2)  J* 


dri2 


(t|— Ti2(t2-n2)(a2-il2) 


Using  equation  (29),  we  can  now  write 
->&2  d^2  -v/ClO2— TJ  )  (r|— EL2  ) 


G(x,t)  )  =  ^ 


(n-’l2^  '/lll('b2-1l2)(a2-Tl2) 


<*> 


Now  using  the  equality 

(Acp) 


=  r  =  i 

T.E.  2 


G(x,T])dT^ 


it  is  possible  to  derive,  after  some  manipulation,  the  relation 


r  =  -tab  ^  i  k»2  K  -  E  +  k^  -  (tf  -  \L  b2k»2  (32) 

where  K  and  E  are  complete  elliptical  integrals  of  the  first  and 
second  kind,  respectively,  of  modulus  k  ±  —  kf 2  = 
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For  the  Kutta  condition  to  apply  at  the  trailing  edge,  equation  (32) 
mu6t  reduce  to  the  expression  r  =  —  yT>E>  p  =  _  ap  in  the  interval 

— t<y<t  (region  II,  fig.  5).  This  leads  to  the  equation  for  the  trail¬ 
ing  edge  in  region  II  given  by 


2kE 

rtk*2 


(33) 


\ 


Figure  7  shows  a  plot  of  the  trailing  edge  derived  from  equation  (33). 
The  asymptotic  slope  of  this  edge  is  45°  and  the  asymptotic  width  of  the 
wing  (i.e.,  b-a  as  b— ^  »)  can  be  shown  to  be  l/n.  By  nfeans  of  this 
plot,  it  is  possible  to  find  the  relation  between  aspect  ratio  and  semi- 
span  s,  for  wings  of  the  type  shown  in  figure  5.  This  was  done  numeri¬ 
cally  and  the  results  are  shown  for  two  different  tip  conditions  in  fig¬ 
ure  8. 


Begion  III  (*^s).—  Consider  for  a  moment  that  the  wing  plan  form 
ends  with  region  II;  that  is,  the  wing  has  the  form  designated  as  type 
b  in  figure  8.  The  downwash  in  any  plane  behind  the  straight  portion 
of  the  trailing  edge  of  such  a  wing  is  determined  by  substituting  equa¬ 
tions  (31)  and  (28b)  into  equation  (23)  and  integrating.  Designating 
this  value  of  wQ  as  vm,  there  results  for  t  <  |  y|  <  s 


and  for  0  <  |y  j  <  t 


Woo  =  3V 


>  (3*0 


v«o  =  3^ 

1  +  ( 

Ei_  i\ 

rtt  J 

/ t2— y2  ] 

_  P 

if 

Ki  Sx  (\ 

t)  -  Ex  Fx  ( 

?)1 

L  \ 

/  s  -y  J 

L  V1 

V  \ 

V  J 

where  .Ex  and  Fx  are  the  incomplete 

modulus  kx  =  t/s  and  argument  y/t,  and  Kx 


elliptic  integrals  with 
is  the  complete  form  of 


is  turning  as  if  the  wing  were  continued.  Hence,  any  surface  behind  the 
straight  part  of  the  trailing  edges  of  the  type  b  wing  will  support 
zero  loading, provided,  of  course,  this  surface  is  also  rotating  at  the 
rate  w  =  py.  Thus  the  solution  for  the  type  a  wing  is  exactly  the 
same  as  that  for  the  type  b  wing  with  the  added  condition  that  the 
loading  is  zero  over  portions  of  the  wing  for  which  x>s. 


Results  for  entire  wing.—  The  results  for  the  entire  wing  plan 
form  can  he  summarized  in  the  form  of  the  total  rolling  moment  and 


Ik 
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subsequently  the  value  of  Cj  .  The  rolling  moment  for  both  types  of 
wings  is  found  from  the  expression 


L  =  2pV0  C2  tan2  9 


(35) 


Using  the  solution  for  T  given  by  operating  on  equation  (31) ,  that  is 
for  0  <  |  y |  <  t 


r  =  -  ny 


and  for  t  <  |y|  <  s 


in  equation  (35)  and  reversing  the  order  of  integration,  the  value  of  the 
moment  can  be  obtained  in  coefficient  form  as 


3czp 
0  tan  9 


I - tl—  (l  rtki'4  +  - iEi-iki'2— ^  (36) 

k  0  tan  0  \  8  6s  3  s  / 


where  again  0  has  been  introduced  as  a  multiplier  of  spanwise  and  verti— 

cal  lengths  and  where  kx  =  t/s  =  Jl  -  kx‘2.  Equation  (36)  is  plotted 
in  figure  9  for  the  plan-form-aspect-ratio  relation  given  in  figure  8. 

When  kx  is  zero,  equation  (36)  reduces  to  equation  (21)  for  the  triangu¬ 
lar  wing. 


Arbitrary  Plan  Form 

If  all  the  terms  in  equation  (28a)  are  considered  then  equation  (28b) 
becomes 

m  n— e 

G*(tii)  =  -  h  ^  aBani  2 
n=l 

and,  using  equation  (29),  equation  (30)  can  be  written 


(28c) 
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-M-eB, 


a2 

'  dTU 

n-»u 


Gi  (^Tn)dm 

n-*u 


Now  if 


HI  rr 

„ ,  \  ^  2ll  T”*  n.— 1  n  &  n— 1 

H(x)  =  2|i  —  ~~  \  nBna  /  sn  u  du 


then 


Gi(x,tj)  =  §-*/(b2 


,2  (^2)  N/(b2-^2)(T)2-a2) 


o2  n-2 

~^nBn  r  Hi2  dni 


and,  finally , 


r  =  ir  Y nBn&1 


(l+k2)snn  1u-2k2snI1+1u  du  —  |i  \  Bna*1  — 


JUi 

2H  f  b2k*2  +  ^jYL  ]T  DBnan~1 


n— 1 

sn  u  du 


In  order  to  satisfy  the  Kutta  condition,  r  in  equation  (39)  must  equal 
m 

-41  ^  Bnan.  Thus,  the  equation  of  the  trailing  edge  becomes 


2k  V  ,  n— iT 
S  =  /.  ****  In 


and,  using  equation  (35)>  the  expression  for  rolling-moment  coefficient 
becomes 


zp  Ar  n— 1  n— 1  /  2,  l-3ki  \  3*^* 

- T -  )  ^n^  s  (kiIn+2  +  -  Tn) - (^l) 

P  tan  0  4sp  tan  6  \  n+  2  J  32 

n=l 


1 6 
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where 


In 


Values  of  the  first  six  ^  are  given  in  Appendix  B. 

Since  the  trailing  edge  must  pass  through  the  point  a  =  0  when 
x  =  1,  the  value  of  Bi  must  always  be  1.  The  other  coefficients, 
however,  are  at  our  disposal.  These  can  be  determined  by  specifying  the 
points  through  which  the  trailing  edge  must  pass  and  setting  up  a  series 
of  simultaneous  equations  using  equation  (^0). 

Two  example s  have  been  considered  which  represent  a  wing  with  very 
little  taper  (fig.  10(a))  and  a  wing  which  tapers  to  a  point  (fig.  10(b)) 
The  trailing  edge  for  the  former  of  these  wings,  type  c,  was  found  by 
finding  the  values  of  B2  in  equation  (ko)  which  made  a  =  0.4  when 
b  =  1.2  and  all  the  B*  s  with  higher  subscripts  were  zero.  The  trailing 
edge  of  the  latter  of  the  two  wings,  type  d,  was  determined  by  finding  B2 
B3,  B4,  and  B5  such  that  the  edge  passed  through  four  points  which  were 
on  a  straight  line  joining  the  root  with  the  tip.  Figure  11  shows  the 
actual  edge  resulting  for  the  type  d  wing  together  with  the  intended 
straight  line  and  the  points  chosen.  The  slight  scallop  indicated  in  the 
trailing  edge  was  found  in  the  case  of  the  triangular  wing  (i.e.,  when 
this  method  was  used  to  solve  the  delta-shaped  wing  and  the  re  suits 
compared  with  the  exact  solution  of  the  linearized  equation) to  give  a 

value  of  - about  5  percent  too  low.  However,  this  was  considered 

p  tan  0 

to  be  sufficiently  accurate  to  establish  the  trend  which  is  shown  for 
this  type  of  wing  in  figure  12. 


Ames  Aeronautical  Laboratory, 

National  Advisory  Committee  for  Aeronautics, 

Moffett  Field,  Calif.,  July  22,  19^9. 

APPENDIX  A 

TABLE  OF  IMPORTANT  SYMB0IB 

radius  of  cylindrical  body  in  transformed  space 

(in  portion  of  report  titled  "Triangular  Wing  with 
Body”  ) 

equation  of  trailing  edge  of  swept  wing  in  transformed 
space  (in  portion  of  report  titled  "Svept-Back  Wing, 
see  equation  (25)) 


a,  a(x) 


a,  a(x) ,  a*(y) 
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IT 


A  (X) 

Ar 

t>,  b(x) 
B(X) 

C 

CZ 


Cl 


P 


radius  of  cylindrical  body  or  equation  of  trailing  edge 

reduced  aspect  ratio,  J  Jl-M02  |  times  aspect  ratio 
equation  of  wing  leading  edge  in  transformed  space 
equation  of  wing  leading  edge 
root  chord  of  wing 


rolling  moment 
(2qS)(wing  area) 


damping-in-roll  coefficient 


c)q  -1 

.  S(pS/VG)  - 


( 


Ei(x) 


G(x,n) 

Ki 


Fi(x) 


k 

ki 

l 

L 

Mo 


damping— in— roll  coefficient  of  wing  alone 


complete  elliptic  integral  of  the  second  hind  with 
modulus  kx 


function  introduced  in  equation  (24) 


complete  elliptic  integral  of  the  first  kind  with  modu¬ 
lus  kx 


incomplete  elliptic  integral 

a 

b 

t 

s 

lift 


dt 


J  (l-k!2t2)(l— t2) 


rolling  moment 

Mach  number  of  free  stream 
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P 

q 

r 

E 

s 

S 

sn(u) 

T, t 

vo 

U,  V,  W 


Ut,  vt,  vt 

X,  y,  Z 
X,  Y,  Z 
P 

AP 

q 

r 

n 

9 

P 

i 

p 

4 


rate  of  roll,  radians  per  second 
dynamic  pressure  ^ApV02^ 

radial,  length  in  polar  coordinates 
ratio  of  tody  diameter  to  wing  span 
wing  semispan  in  transformed  space 
wing  semispan 
.Jacotian  elliptic  function 

distance  used  in  study  of  swept-tack  wing  (see  fig.  5) 
free-stream  velocity 

transformed  perturbation-velocity  components  in  x,y,z 
directions,  respectively 

radial  component  of  velocity 

perturtation— velocity  components  in  X,Y,Z  directions, 
respectively 

transformed  Cartesian  coordinates 
Cartesian  coordinates  in  physical  plane 

J\  l-«o2| 

loading  coefficient, 

/■pressure  on  lower  surface  -  pressure  on  upper  surface 

V  q 

spanwise  distribution  of  circulation 

y2 

semiapex  angle  of  wing 
pC2  tan2  9 

complex  variable  (y  +  iz) 
free— stream  density 
perturbation  potential 
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perturbation  potential  in  transformed  coordinates 
jump  in  potential  in  the  z  =  0  plane 
polar  angle 

APPENDIX  B 

Exact  expressions  for  Ij.,  I2,  .  .  .  I6* 


9 

AcP0 

t 


Ii  =  E 


_  1  k'2  7  1+k 

I2  =  —  +  -  In  - 

24k  1-k 


Ia  =  3^-  [  k»2K  -  (l-2k2)  E] 


T  =  1  _  1  k^f_ 

4  8  k2 


1  _  m 

2k  1-k 


) 


I*  = 


^  2(l-4kg)(l+kg) 


15k2 


K-E  1— 4k2 

k2  15k2 


K 


(l+3k2)(3-5k2)  ( 1-k2 ) ( l+2k2+5k4 )  1+k 

j  = - - -  +  - r - 


6 


48k' 


32k 


1-k 


Notice  that  when  k=l,  In  = 


n 


Series  expansions  for  Ii,  l2> 

Ii 


-i(l-£ 

2  \  4 


.  I6  for  k  <1. 


3  ^4 - k&_ 


64 


256 


...) 


Ip  =  1  -  -  k2  -  —  k4  -  —  k6- 

2  3  15  35 


j°  =  l  C1-!*2 


5  k4 - ^-r  k6 


64 


1024 


2  __ 

i-k2- 

-i-k4- 

-1  ka.  .  .  . 

3 

15 

35 

315 

rt 

(3 

k2--^ 

k+ - —  k6—  .  . 

16 

\  4 

128 

512 

_8 

-  _§  k2 

-  ii-  k4 

16  k6 

15 

35 

315 

693 
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Figure  I.-  Comparison  of  two  theoretical  results  for 
rolling  triangular  wing. 
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Figure  4-  Effect  on  damping  in  roll  of  a  cylindrical 
body  mounted  on  a  pointed  wing  with  unswept 
trailing  edge. 


Transformations 
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swept- back  wing. 


NACA  TN  1950 


27 


of  trailing  edge  y-  a(x)  for  special  swept- back  wing. 
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for  low  values  of  reduced  aspect  ratio. 


IS  -  Variation  of  damping-in-roll  parameter  with  reduced 
various  wing  plan  forms  and  for  low  values  of  reduced  < 


